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ABSTRACT 

Whilst young massive clusters (YMCs; M > 10"^ Mq, age < 100 Myr) have been 
identihed in significant numbers, their progenitor gas clouds have eluded detection. 
Recently, four extreme molecular clouds residing within 200 pc of the Galactic centre 
have been identified as having the properties thought necessary to form YMCs. Here 
we utilise far-IR continuum data from the Herschel Infrared Galactic Plane Survey 
(HiGAL) and millimetre spectral line data from the Millimetre Astronomy Legacy 
Team 90 GHz Survey (MALT90) to determine their global physical and kinematic 
structure. We derive their masses, dust temperatures and radii and use virial analysis 
to conclude that they are all likely gravitationally bound - confirming that they are 
likely YMC progenitors. We then compare the density profiles of these clouds to those 
of the gas and stellar components of the Sagittarius B2 Main and North proto-clusters 
and the stellar distribution of the Arches YMC. We find that even in these clouds - 
the most massive and dense quiescent clouds in the Galaxy - the gas is not compact 
enough to form an Arches-like {M = 2x10'^ Mq, Re// = 0.4 pc) stellar distribution. 
Further dynamical processes would be required to condense the resultant population, 
indicating that the mass becomes more centrally concentrated as the (proto)-cluster 
evolves. These results suggest that YMC formation may proceed hierarchically rather 
than through monolithic collapse. 

Key words: Stars: formation - ISM: clouds - Galaxy: centre, open clusters and 
associations: general 


1 INTRODUCTION 


Stellar clusters can contribute substantially to the stel¬ 


lar masses and densities can reach and even exceed those of 


lar population of a galaxy (e.g. 

Lada & Lada 2003 [God- 

dard, Bastian & Kennicutt|2010| 

Silva-Villa & Larsen|2011| 

Kruijssen|2012 

1 . Classically, they have been categorised into 


two distinct populations; globular clusters (GCs) and open 
clusters (OCs). More recently however, it has become appar¬ 


ent that young massive clusters (YMCs; defined by Portegies 
|Zwart, McMillan fc Gieles|2010| to be gravitationally bound 
systems with M > 10“^ M© and ages < 100 Myr), whose stel- 


E-mail: D.L.Walker@2009.1jmu. ac.uk 


GCs, are still forming at the present day (e.g. Holtzman 
|et al.|1992[ ). This discovery has reinvigorated the discussion 
of cluster formation mechanisms and spawned the idea that 
perhaps all clusters form in a similar way, with OCs rep¬ 
resenting the low-mass counterparts and GCs representing 


the early-Universe analogues of extreme YMCs (Elmegreen 
fc Efremov||1997 Kruijssen]2014 1. Understanding and char¬ 


acterising the formation and evolution of YMCs may there¬ 
fore be fundamental in revealing how clusters form across 
the full mass-range. 

YMCs with ages greater than a few Myrs have been 
identified in significant numbers (see [Portegies Zwar^ 
McMillan fc Gieles|2010 and references therein), and partic- 
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Figure 1. Top: Masked column density map of the ‘dust-ridge’, derived from HiGAL far-IR continuum observations ( [Battersby et al~] 
|201l[ |. A column density threshold of 5 x 10^^ cm“^ (red contours) is implemented and contours (white) are at 10, 15 ... 55 x 10^^ cm“^. 
Bottom: Masked temperature map of the ‘dust-ridge’, derived from HiGAL far-IR continuum observations (for the techniques used to 
derive these maps, see [Battersby et al.||20fT| . Contours (white) are at 17, 18 ... 27 K (systematic uncertainties in SED fitting are ~ 
3—4 K). To display the correlation between dust column density and temperature, the same pixels as in the top image are masked here. 
Clearly seen in all of the clumps is an anti-correlation between column density and temperature of the dust — consistent with them being 
centrally condensed with cold interiors ||Longmore et al.|2012^. 


ularly in galaxies with high star-formation rates (e.g. Whit- 
more|2002 1. Several have also been identified within our own 
Galaxy, such as the Arches, Quintuplet, NGC 3603, Wester- 
lund 1 and red super-giant (RSG) clusters (Figer et al.|1999 


Goss fc Radhakrishnan|1969[[Clark et al.|2005| [Davies et al.| 


2007). However, if we are to understand how such clusters 


form, we must identify and characterise their progenitor gas 
clouds. These progenitor clouds are expected to contain > 


10“ Mq within only a few parsecs. Indeed, Ginsburg et al. 


(20121 and Urquhart et al. (20131 find several such clouds 


in the Galactic disk that satisfy these conditions. However, 
these clouds are all actively forming stars and will there¬ 
fore have lost much of their initial structure as a result of 


stellar feedback (note that Urquhart et al. (20131 explic¬ 
itly searched for star-forming clouds). Whilst these clouds 
are of course important in understanding the YMC forma¬ 
tion process, if we are to probe the initial conditions then 


we require clouds that have yet to initiate widespread star 
formation. Thus, in addition to the aforementioned crite¬ 
ria, we would then also expect such progenitor clouds to be 
cold (T ~ 10 - 20 K) and devoid of widespread, high-mass 
star formation. Given that none of the identified progenitor 
clouds in the Galactic disk are quiescent, very little is known 
about the initial conditions from which these massive clus¬ 
ters form. Identifying and characterising the gas within such 
clouds is therefore vital if we are to understand the process 
of YMC formation. They also present a unique insight into 
the earliest stages of star, planet and cluster formation in 
extremely dense environments, where feedback mechanisms 
and dynamical encounters are likely to affect the subsequent 
cluster evolution significantly. 


Longmore et al. (20121 recently identified a likely 


YMC precursor in the extreme infrared dark cloud (IRDC) 
GO.253-1-0.016 (see e.g. Lis et aL]|1994 Lis fc Menten|[r998| 
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Longmore et al.||2012| 

Kauffmann, Filial Zhang] 

20131 

Rathborne et al. [[2014a 

Johnston et al.[[2014[ Rath 

3orne[ 


and determined that it is cold (T^ust ~ 20 - 
30 K), high-mass (M ~ 1.3 x 10® M©), compact (r ~ 2.8 pc) 
and other than a single region of weak water-maser emission, 
exhibits little-to-no star forming activity - precisely the con¬ 
ditions one would expect of a dense cloud that could form 


a high-mass cluster. Longmore et al. (2013b I later identified 


three further potential YMC precursors in clouds known as 
‘d’, ‘e’ and ‘f’. Immer et al. (20121 show that these clouds 

pc-scale) and 


are all high-mass (~ 10® Mq), compact ( 
other than a region of methanol maser emission towards 
cloud ‘e’, are all quiescent. Along with GO.253-1-0.016, they 
are situated in the Central Molecular Zone (CMZ; inner ~ 
200 pc of the Galaxy, Morris fc Serab>^ 1996[) and belong 


to the so-called ‘dust-ridge’ (Lis et al. 1999[ [Immer et ah 


derive these maps). A column density threshold has been ap¬ 
plied such that all pixels below 5 x 10^^ cm“^ are masked. 
This threshold is chosen as it highlights the spatial extent of 
the dense clumps well, whilst effectively masking the more 
diffuse emission across the region. GO.253-1-0.016, ‘d’, ‘e’ and 
‘f’ have peak column densities ranging from 2.6 - 5.3 x 10^® 
cm“^ and central temperatures ranging 17 - 19 K. Over¬ 
all, the maps clearly display an anti-correlation between the 
column density and temperature of the dust in every cloud 
- consistent with them being centrally condensed with cold 
interiors. We note that the central temperature of cloud ‘c’ 
is slightly higher than might be expected given its central 
column density - at 22 K it is ~ 4 K warmer in its core than 
the other clouds, despite all being at similar column densi¬ 
ties. This is likely due to the fact that this cloud is likely 
forming high-mass stars that heat the cloud’s interior, as 


2012) towards the Galactic centre, which itself appears to evidenced by water and methanol maser emission (Forster 


belong to a coherent circumnuclear gas ring. |MoIinari et al.| 
(2011) model this ring as a twisted, oo-shaped ellipse that 


orbits the central super-massive black hole, Sgr A*, with 
a semi-major orbital axis of ~ 100 pc. [Kmijssen, Dale~^ 


Longmore (20151 confirm the coherence of the gas stream 


using dynamical modelling. The fact that four of the most 
extreme clouds known to exist in the Galaxy all lie at the 
same distance and reside within the same stream of gas is 
certainly very intriguing, but it also presents an opportunity 
to study a sample of potential YMC progenitor clouds un¬ 
der the same environmental conditions at similar sensitivity 
and resolution. 

In this paper we extend the analysis of |Longmore et J?] 
( |2012[ ) (L12, hereafter) of GO.253-1-0.016 to include clouds 
‘d’, ‘e’ and ‘f’. We utilise continuum data from the Her- 
schel infrared Galactic Plane Survey (HiGAL,|Molinari et al. 


[Mo] ) to measure their global physical properties such as 
mass, radius and temperature. These HiGAL data cover a 
wavelength range of 70 - 500 /rm and provide ~ 5.5” or 
0.22 pc resolution at the distance of these clouds, which is 
taken to be ~ 8.4 kpc ( Reid et al.||2009 ). Spectral line data 
from the Millimetre Astronomy Legacy Team 90 GHz Sur¬ 
vey (MALT90, [Foster et al.|2011[|201^ [Jackson et al.|2013[ ) 
are used to measure the global kinematic properties of these 
clouds. The MALT90 survey provides spectral line data for 
16 lines in the 90 GHz band. The data cubes for these tran¬ 
sitions consist of 4096 channels with 0.11 km s“^ velocity 
resolution. Using the results from these data, we determine 
whether these clouds - the most massive and dense quies¬ 
cent clouds known in our Galaxy - are capable of forming 
YMCs. Having shown that they are indeed candidate YMC 
precursors, we then compare these clouds to more evolved 
(proto)-YMCs to speculate as to how gas is converted to 
stars on large scales in the early stages of YMC formation. 


2 RESULTS 

2.1 Dust Column Densities &: Temperatures 

Figure 1 {upper panel) shows the HiGAL column density 
map of the ‘dust-ridge’ at the Galactic centre. The bottom 
panel displays the HiGAL temperature map of the same re¬ 
gion. (See Battersby et al.[[2011 for the techniques used to 


& Caswell 1999 Caswell et al. 2010). It is interesting to 


note the discrepancy between the gas and dust tempera¬ 


tures at the Galactic centre. Ao et al. (20131 find that the 


gas temperature towards the Galactic centre ranges from 50 
K to in excess of 100 K. This is significantly higher than 
the low dust temperatures of ~ 20 K, suggesting that the 
gas is being heated by some non-photon driven mechanism 
such as cosmic ray heating or turbulent energy dissipation. 
With respect to our sample of clouds, we note that [Ao et al.[ 
120131 only measured the temperature of GO.253-1-0.016. We 


assume that their results also hold for the other three clouds. 


2.2 Dust Mass 


Using these column density maps, we estimate the dust 
masses of the clouds by assuming a mean molecular weight of 
2.8m_ff, multiplying each pixel by its physical area (assuming 
a distance of 8.4 kpc) and summing over the cloud area. For 
GO.253-1-0.016 and cloud ‘d’, this is straightforward given the 
column density threshold of 5 x 10^^ cm“^. However, clouds 
‘e’ and ‘f’ are not entirely distinct, are embedded in a higher 
density region that lies above this threshold, and have two 
nearby distinct clumps at different velocities. Any kinemat¬ 
ically distinct emission (identified using MALT90 data; see 
^ 2 . 4 ) is therefore masked to determine the masses of clouds 
‘e’ and ‘f’. We determine the masses of GO.253-1-0.016, ‘d’, 
‘e’ and ‘f’ to be 11.9 x 10^ Mq, 7.6 x 10^ Mq, 11.2 x 10"* 
Mq and 7.3 x lO'* Mq, respectively (see Table 1). We note 
that these mass estimates are necessarily uncertain due to 
the ambiguity involved in defining a distinct cloud area, par¬ 
ticularly in a contiguous region such as this. 

Figure 2 shows how the estimated masses decrease as 
we choose higher column density thresholds. Increasing the 
threshold from 5 - 10 x 10^^ cm“^ decreases our mass es¬ 
timates by approximately 20% for GO.253-1-0.016, 40% for 
cloud ‘d’, 10% for cloud ‘e’ and 15% for cloud ‘f’. Given 
the systematic uncertainties in estimating dense gas mass 
towards the Galactic centre of a factor of ~ 2 (Longmore 


et al. 2013al, any uncertainty in these mass estimates is 


dominated by systematics rather than the threshold used. 
These estimates agree well with those found by [Immer et al.[ 
(2012) and L12. Note that the slight (<10%) discrepancy 
between L12’s and our mass estimate for GO.253-1-0.016 is 
simply due to a slightly higher column density threshold. 
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Figure 2. Mass determination as a function of column density 
threshold for clouds GO.253+0.016 {green), ‘d’ (black), ‘e’ (red) 
and ‘f’ (blue). We implement a lower threshold of 5 x 10^^ cm“^ 
and incrementally increase this to 10 x 10^^ cm“^. Clouds ‘e’, ‘f’ 
and GO.253+0.016 are not overly sensitive to the column density 
threshold, with ~ 10% and ~ 15% changes over this range, respec¬ 
tively. The mass of cloud ‘d’ is however less well constrained, with 
a change in mass of ~ 40% over the range. The legend displays 
the mass determined for each cloud at the lowest threshold. 

2.3 Radius 

Defining characteristic cloud radii is further compli¬ 
cated by their non-spherical geometry. Rather than simply 
prescribing representative circular radii, we instead take a 
geometric mean of the minor and major axes of each cloud, 
whose boundaries are defined by the aforementioned column 
density threshold, to derive effective radii. We derive effec¬ 
tive radii of 2.9 pc, 3.2 pc, 2.4 pc and 2.0 pc, respectively, 
for GO.253+0.016, ‘d’, ‘e’ and ‘f’. 

To compare how compact these clouds are, we deter¬ 
mine characteristic radii within which 5 x 10"^ Mq is en¬ 
closed. This mass is chosen as we know that all four clouds 
exceed this mass and it is large enough such that a 10 "^ Mq 
cluster could form from it, given a star formation efficiency 
as low as 20%. Figure 3 displays this mass-radius relation for 
the four clouds and black dashed lines indicate their charac¬ 
teristic radii - 1.0 pc, 1.2 pc, 0.8 pc and 1.0 pc, respectively 
for GO.253+0.016, ‘d’, ‘e’ and ‘f’. This highlights the extreme 
nature of these clouds in that they harbour enough mass to 
form a YMC within a radius of ~ 1 pc and yet they do not 
show any signs of high-mass star formation other than that 
inferred by the detection of weak maser emission. 

2.4 Gas Velocity Dispersion 

We use molecular line emission to determine the global 



Figure 3. Enclosed mass as a function of radius for GO.253+0.016 
(green), ‘d’ (black), ‘e’ (red) and ‘f’ (blue). Intersections of black 
dashed lines indicate the radii within which 5 x 10^ Mq are en¬ 
closed (see legend). 

kinematics of these clouds and assess their dynamical state. 
This analysis is complicated due to their location in the 
CMZ, where the kinematic and chemical structure is known 
to be complex - with large velocity dispersions, shock- 
enhanced chemistry and confusion due to unrelated line- 
of-sight emission within the Galactic disk. It is therefore 
important to ensure that we use sensible molecular tracers 
to isolate the emission from each cloud. 

The MALT90 survey provides data cubes for 15 differ¬ 
ent molecular tracers, probing a range of critical densities 
and excitation energies, thus making it ideal for our pur¬ 
poses. Across clouds ‘d’, ‘e’ and ‘f’, emission from 11 of 
the 15 tracers is well detected. We adopt the same cat¬ 
egorisation as [Rathborne et ah] ( |2014a[ ), separating these 
detected molecular transitions as: optically thick tracers - 
HGN (1-0), HNC (1-0), HGO+ (1-0) and N 2 H+ (1-0); op¬ 
tically thin tracers - ( 1 - 0 ), H13G0+ ( 1 - 0 ) and 

(2-1) and ‘hot core’/shock tracers - HNGO 4(0,4)-3(0,3), 
SiO (1-0), HC 3 N (10-9) and GH 3 CN 5(0)-4(0). As might 
be expected of clouds in the turbulent environment of the 
CMZ, we find that their kinematics are complicated. The 
data cubes show complex line profiles, significant velocity 
gradients, multiple velocity components, large line-widths 
and intense shock-associated emission lines. We defer de¬ 
tailed analysis of these features to a subsequent paper. (See 
[Rathborne et al)] ( |2014a| ) for an in-depth discussion regard¬ 
ing the molecular line emission from GO.253+0.016). Here 
we simply wish to isolate the emission from each cloud to 
obtain estimates of the line-widths. To do this, we make 
the reasonable and common assumption that the optically 
thin transitions trace the underlying kinematics of the clouds 
most reliably, whereas optically thick lines only trace the 
surface kinematics of the clouds. Of these, HN^^C is the 
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brightest line detected in all three clouds and is therefore 
used for line-width determination in all cases. Figure 4 shows 
Hanning-smoothed HN'^^C spectra that have been averaged 
over the spatial extent of the clouds (defined by the 5 x 10^^ 
cm“^ column density threshold, see Figure 1), where each 


profile is fitted using the CASA software package (|McMuUTn 


et al.|2007 |. Cloud ‘d’ {left panel) displays two velocity com¬ 
ponents, one at ~ 0 km s“^ and another at ~ 20 km s“^. 
We attribute only the 20 km s”'^ component to the gas in 
cloud ‘d’, as the morphology of the molecular line emission 


ponent arises from an unrelated cloud along the line of sight 
(see Figure 5a). The fit to the associated component yields 
a peak velocity and line-width (FWHM) of V lsr = 20.75 ± 
0.69 km s“^ and AV = 16.3 ± 1.5 km s“^. Cloud ‘e’ {middle 
panel) shows a clear singular component with Vlsr = 31.66 
± 0.22 km s“^ and AV = 15.51 ± 0.52 km s“^. Cloud ‘f’ 
{right panel) also shows two velocity components, one at ~ 
27 km s“^ and another at ~ 40 km s“^. Though both com¬ 
ponents show similar line-widths and peak intensities, we 
conclude that only the 40 km s“^ component is associated 
with cloud ‘f’, as again the molecular line emission matches 
the dust emission well (see Figure 5b). This component has 
a peak velocity and line-width of Vlsr = 40.1 ± 2.6 km 
s“^ and AV = 16.5 ± 3.2 km s“^. Coupling these results 
with L12’s results for GO.253-1-0.016 - Vlsr = 36.1 ± 0.4 
km s“^, AV = 15.1 ± 1.0 km s“^ - we see that all four 
clouds, whilst having slightly different peak velocities, have 
very similar line-widths of ~ 16 km s“^. 


2.5 Virial mass — Gravitationally bound? 


If we are to investigate the cluster-forming potential of 
these clouds, we must first assess whether they are gravita¬ 
tionally bound. To do this we estimate and compare their 
dust and virial masses. Taking the radii and line-width es¬ 
timates and assuming a uniform spherical density distribu¬ 
tion, we use M„ir = fcRAV^ (where k = 126 and the units 
of M, R and AV are Mq, pc and km s“^, respectively; see 


MacLaren, Richardson & Wolfendale 19881 to derive virial 
masses of 1.1 x 10° M©, 7.3 x 10'‘ M© and 6.9 x 10"^ Mq for 
clouds ‘d’, ‘e’ and ‘f’, respectively. 

Under these assumptions, we compare the dust masses 
with the virial masses to obtain virial ratios of a = 
Mvir/-A4dust = {1.4,0.7,0.9}. It is clear that to within ~ 
50%, Mvir ~ Afdust in all cases. We therefore conclude that 
all three clouds are likely gravitationally bound. We note 
that the effect of magnetic pressure has not been accounted 
for in this analysis. 



Figure 5. (a) Integrated HNCO intensity maps of cloud ‘d’ show¬ 
ing emission integrated between -10 to 5 km s“^ {left) and 5 to 50 
km s“^ {right), (b) Integrated HNCO intensity maps over clouds 
‘e’ and ‘f’ showing emission integrated between 5-25 km s“^ 
{left) and 25 - 50 km s“^ {right). Overlaid in white are contours 
from the HiGAL column density maps. These maps show that 
the emission over the lower velocity range does not match well 
the morphology of the dust emission toward these clouds. Thus, 
we assume that this emission arises from unrelated clouds along 
the line of sight. Note that HNCO emission is used here as it 
matches the HN^^C emission very closely but is detected with a 
much higher signal to noise. 


has only been measured for GO.253-1-0.016 - we are assuming 
that this is also true for the other clouds. Any associated un¬ 
certainty in the sound speed is smaller as Cs ~ T®'®). This 
is then used to estimate the Mach number as A4 = cr/cs 
(where cr = AU/2.355). We estimate the sound crossing time 
as tsc = R/cs and the cloud crossing time as tec = R/c. Fi¬ 
nally, we estimate the free-fall time as ts = ^JTTJYTGp, 
where p has been estimated under the assumption of a uni¬ 
form spherical distribution. 


3 DISCUSSION 


2.6 Time-scales 

Using the derived masses, radii and line-widths, we de¬ 
rive the sound crossing times, cloud crossing times and free- 
fall times for these clouds (see Table 1). Firstly we estimate 
the sound speed as Cs = -y/fcTg/p, where k is the Boltz¬ 
mann constant, p is the mean molecular weight (taken to 
be 2.8mH) and Tg is the gas temperature (assumed to be 70 
K in all cases, [Ao et al.|r2013[ though we reiterate that this 


Table 1 displays the derived global properties for clouds 
GO.253-1-0.016, ‘d’, ‘e’ and ‘f’, along with characteristic den¬ 
sities and time-scales for each cloud. These results highlight 
that all of these clouds are extremely massive and compact 
objects, containing ~ 10 ® Mq within only a few parsecs, 
resulting in correspondingly high column and volume densi¬ 
ties and short dynamical time-scales. Goupling these results 
with the fact that the clouds are close to virial equilibrium, 
the conclusion of |Longmore et ^ ( |2013b| ) is confirmed; that 
these clouds are all excellent YMC progenitor candidates. 
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Figure 4. Hanning-smoothed, spatially-averaged HN^^C line profiles for clouds ‘d’ [left], ‘e’ [centre] and ‘f’ [right]. Each is fitted using 
the multi-component Gaussian-fitting in the CASA software package. In the profile for cloud ‘d’, we attribute the component at ~ 0 km 
s~^ to emission from an unrelated cloud along the line of sight and in cloud ‘f’ we attribute the ~ 27 km s~^ component to emission 
from nearby clouds (see Figure 5 for integrated intensity images). Omitting these unassociated components, we obtain the following peak 
velocities and line-widths: cloud d: Vlsr = 20.75 di 0.69 km s“^, AV = 16.3 di 1.5 km s“^; cloud e: Vlsr = 31.66 di 0.22 km s“^, 
AV = 15.51 ± 0.52 km s-^; cloud f: Ylsr = 40.1 ± 2.6 km s"!, AV = 16.5 ± 3.2 km s-A 


Intriguingly, there are no clouds yet known to exist in the 
rest of the Galaxy that have such extreme masses and densi¬ 
ties that are not forming stars prodigiously ( [Ginsburg et al.| 


forms in an extremely compact natal gas cloud (i.e. 
Ri^it < Rfin, plZt > p/in) - Feedback processes grad¬ 
ually remove the remaining gas, diluting the global 
potential and causing the cluster to expand towards 


and dense as these currently exhibit minimal signs of active 

its final, un-embedded phase (see e.g. 

Lada, Margulis 

star formation may therefore have important implications. 

^ Dearborn||1984| |Boily & Kroupa||2003[ |Bastian ^ 

For example, it may suggest that star and cluster formation 

Goodwin|2006 

Baumgardt & Kroupa|2007 

1 . This re- 


proceeds differently at the Galactic centre than in the disk 
(see |Longmore et al.||20f4| and references therein), or that 
we have simply caught the earliest stages of YMC forma¬ 
tion in four separate cases. Given that the environment at 
the Galactic centre is extreme compared to the disk ( |Shetty| 
et al.|2012 Kruijssen fc Longmore|2013 1, it is plausible that 


star formation may be inhibited in some way (Longmore 
|et al.|[2013a| [Kmijssen et al.||2014[ ). However, it is known 
that YMCs can and do form near the Galactic centre since 
at least two YMGs, the Arches and Quintuplet clusters, have 
formed there. The existence of YMCs at the Galactic centre 
therefore adds further weight to the conclusion that these 
extreme clouds have the capacity to form YMGs. We now 
investigate how these clouds can be used to gain insight in 
to the process of YMC formation. 


suits in a bound, spherical cluster with Rfi„ > 

(ii) Stars and sub-clusters form in a gas cloud with 
> Rfin, pZit < P/in- They form throughout 
the spatial extent of their natal gas clouds, following 
the hierarchical structure of the interstellar medium 
( Larson||1981 \. A heightened star formation efficiency 
(SFE) in the densest peaks leads to gas exhaustion 
on local scales, causing stellar dynamics to eventu¬ 
ally dominate ( Kruijssen et al.|2012 Girichidis et al. 
|2012[ ). The subsequent hierarchical merging of these 
condensations results in a centrally concentrated, 
bound cluster (Fujii, Saitoh fc Portegies Zwart|2012 


[Parker et al.|[^14[). 


3.1 Comparing Clouds &: Clusters 


Having derived the global properties of GO.253-1-0.016, ‘d’, 
‘e’ and ‘f’ and confirming that they are sufficiently mas¬ 
sive and dense to potentially form YMCs, we now compare 
their properties to the observed properties of the intermedi¬ 
ate and final stages of YMC evolution. The reasoning here is 
that through comparing YMCs in their initial, intermediate, 
and final stages, we can begin to build up a coherent picture 
of how these different stages connect to one-another and ul¬ 
timately determine whether the process of YMC formation 
is fundamentally different from that which forms low-mass 
clusters. We posit two cluster formation scenarios, which can 
be distinguished by comparing the initial gas sizes and den¬ 
sities plnlt) with those of the resultant stellar popu¬ 


lation {R*fi„, p*fin) (see Longmore et al.|2014 and references 
therein for a more detailed discussion): 


(i) A bound, centrally-condensed stellar population 


In essence, we should be able to distinguish between 
these two scenarios by studying the most likely progenitor 
systems - massive and compact molecular clouds. If sce¬ 
nario (i) is a common mode of YMC formation, we ought 
to see > 10® Mq clouds that are more or equally as com¬ 
pact as Galactic YMCs. If instead scenario (ii) is favoured, 
then we might expect to see such clouds with sizes larger 
than those of YMCs that show fragmented sub-structure on 
small spatial scales. Given that GO.253-1-0.016, ‘d’, ‘e’ and 
‘f’ represent the most extreme quiescent molecular clouds 
known in the Galaxy, they offer an ideal sample of progen¬ 
itor systems with which we can investigate the validity of 
these two scenarios prior to the loss of initial structure due 
to feedback from high-mass star formation. We can begin 
to distinguish between these scenarios simply by looking at 
their mass surface density profiles to see how they compare 
to Galactic YMCs and proto-YMCs. 

We use the HiGAL column density maps to obtain mass 
surface density profiles for clouds ‘d’, ‘e’ and ‘f’ by calculat¬ 
ing the enclosed mass (see ^2.S) within increasing circular 
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apertures (centred on column density peaks) and dividing by 
the area of the corresponding aperture. For GO.253+0.016, 


appendix A) - 


we use ALMA Cycle 0 + single-dish data (Rathborne et al. 


2014b I. These data consist of ALMA observations at 3 mm, 
with an angular resolution of 1.7” (0.07 pc at a distance 
of 8.4 kpc). The data were then combined with Herschel 
500 pm to recover the flux that is filtered out by the inter¬ 
ferometric observation. The resolution of the ALMA data 
allows us to investigate the surface density prohle of this 
cloud down to much smaller radii. Interestingly, |Rathborne| 


et al. (2014a I do find evidence for fragmentation and hierar¬ 


chical small-scale structure in GO.253+0.016. Furthermore, 


the ALMA data (Rathborne et al. 2014bI directly reveals 


such sub-structure within this cloud, showing that the gas 
is highly fragmented and contains a number of dense cores 
distributed throughout the cloud (we direct the reader to 
[Rathborne et aL]|2015| for detailed analysis of the cloud’s 
sub-structure). This suggests that the internal structure of 
GO.253+0.016 is consistent with that predicted by scenario 
(ii). Higher resolution data are required for the other clouds 
to assess whether this is true in all cases. 

In all following plots, we have applied a multiplicative 
factor of 1/3 to the data for GO.253+0.016, ‘c’, ‘d’, ‘e’ and 
‘f’ such that we can infer what the resultant stellar popula¬ 
tion would look like if it were to form with a star formation 
efficiency (SFE, e) of 1/3 at the current mass distribution 
of the clouds. In reality, SFE will of course vary throughout 
the cloud, where it will be enhanced towards local density 
peaks - this factor is chosen to represent a global SFE. 


3.1.1 Identifying the Intermediate Phase 

To investigate the active star-forming phase of YMC forma¬ 
tion, we choose to study the gas and stellar content of the 
well-known star-forming Sagittarius B2 complex (Sgr B2, 
e.g. Qin et al. 2011| . We include both Sgr B2 Main and 
North, as these are potential proto-YMCs in a deeply em¬ 
bedded phase and lie in the same region as GO.253+0.016, 
‘d’, ‘e’ and ‘f’. Qin et al. (20111 use data from the Submil¬ 


limeter Array (SMA) to investigate the dense sub-structure 
within Sgr B2 Main and North. Interestingly, they find that 
Sgr B2 Main contains many sub-mm sources that appear 
to have a fragmented spatial distribution, whereas Sgr B2 
North only contains two sub-mm sources. They propose that 
this may suggest that the main cluster is more evolved and 
that northern cluster is less evolved and characterised by 
monolithic high-mass star formation. 

The emission from the gas and dust in Sgr B2 is satu¬ 
rated in the HiGAL data. To obtain surface density profiles 
for the gas in Sgr B2 Main & North, we instead use data from 
the Bolocam Galactic Plane Survey (BGPS, |Rosolowsky| 


et al.|2010 Aguirre et al.|20TT Ginsburg et al.|2013 (. These 

data are at a wavelength of 1.1 mm and provide a pixel scale 
of 7.2” (~ 0.3 pc at a distance of 8.4 kpc). 

As we are using these data to generate mass surface 
density profiles for the clouds, we convert the data from units 
of intensity to units of mass. We do this using the following 
relation (taken from Kauffmann et al.|2008 equation A.31, 


M = O.12M0 




439(A/mm)“^(T/10 K)" 


d 


-5 

mm J 


( 1 ) 


o.oicm^g-iy VJyy Vioopc^ 
where M is mass, A is wavelength, T is the dust tempera¬ 
ture, ku is the dust opacity, 7+ is the integrated flux and d 
is distance. The dust temperature is assumed to be 20 K, 
though we note that a star-forming complex like Sgr B2 will 
not be isothermal due to heating from the embedded HII re¬ 
gions and stellar population. Additionally, it is known that 
the kinematic structure of Sgr B2 is complex, with multiple 
velocity components towards the region. As a result, any 
mass and density estimates for the gas in this region will be 
upper limits. Schmiedeke et al. (in prep) model the region 
in more detail. 

The only observationally unconstrained parameter in 
the above relation is the dust opacity (k,^). To estimate this, 
we use the following relation, given in §3.2 of|Battersby et al.| 

( MT| ) - 

1.75 


k,y = 0.04 cm g 


V505 GHz/ 


( 2 ) 


.505 GHz/ 

where v is the frequency. Note that this contains the explicit 
assumption that the gas-to-dust ratio is 100. 

Obtaining a mass surface density profile of the proto- 
cluster(s) within Sgr B2 is difficult due to the high col¬ 
umn densities and hence extinction towards the region. To 
overcome this, we hrst take the positions and zero age 
main-sequence (ZAMS) spectral classifications of the stel¬ 
lar sources embedded in the Ultra-Compact HII (UCHH) 
regions from Tables 2 and 3 in Gaume et al. (19951. We 


then convert the spectral type of each source to a repre¬ 
sentative mass using the spectroscopic masses of ZAMS OB 


stars given in Table 5 (column 8) of Vacca, Garmany & 


Shull (19961. Knowing the spatial distribution and masses 


of the OB stars within Sgr B2 and assuming a distance of 
8.4 kpc, we then proceed as previously and calculate the 
total stellar mass enclosed within increasing circular aper¬ 
tures. Given that the observations in [Gaume et al.| ( [1995[ ) 
are sensitive only to stars > 10 Mq, we correct for the to¬ 
tal mass by applying Kroupa-type IMF ( Kroupa|2001 1 , for 
which the fraction of mass >10 Mq is ~ 0.16. The total mass 
is therefore estimated by applying a multiplicative factor of 
1/0.16 and normalising by the mean stellar mass of ~ 0.5 
Mq. Taking all 25 sources in [Gaume et al.[ ( [1995[ ), we esti¬ 
mate a total stellar mass of ~ 3.5 x 10'” Mq. The UCHH 
regions in Sgr B2 also indicate a sub-structured distribution 
of the embedded high mass stars - possibly consistent with 
scenario (ii). 


3.1.2 Comparison to Galactic YMCs 

For comparison with Galactic YMCs, we choose to discuss 
the Arches cluster as it is massive (M ~ 2 x lO’^ M©), com¬ 
pact {Reff ~ 0.4 pc), young (Age ~ 2 Myr) and is situated 
towards the Galactic centre and therefore a similar environ¬ 
ment to Sgr B2, GO.253+0.016, ‘d’, ‘e’ and ‘f’ ( [Portegies] 
Zwart, McMillan fc Gieles||2010 K If YMC radii and central 
densities are related to the tidal radius, such that more com¬ 
pact YMCs are formed in stronger tidal fields, it is impor¬ 
tant that we compare clouds and clusters within the Galactic 
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centre so as to eliminate any environmental variations. The 
Quintuplet cluster is also situated towards the Galactic cen¬ 
tre. However, it is an older system than the Arches and it has 
been suggested that the disruption time-scale of clusters at 
the Galactic centre is short, occurring over ~10 Myr (Kim, 


Morris fc Lee|[T999| [Portegies Zwart et al.|[2002[ |Kruijsse~ 
et aI.|20T4 |. As such, we use only the Arches for comparison 


here as it is more likely representative of an initial YMC 
distribution. We utilise the observed surface density profile 
shown in Figure 16 of Espinoza, Selman & Melnick (20091 


along with the given best-fit King density profile and param¬ 
eters, correcting for the number and masses of stars below 
10 Mq with a Kroupa-type IMF. We also use the observed 
cumulative mass profile given in Figure 3 of |Harfst, Porte-] 
gies Zwart & Stolte (20101, again IMF-corrected for stars 
below 10 Mq. 


Figure 6 displays the enclosed mass as a function 
of radius for GO.253-1-0.016, ‘d’, ‘e’ and ‘f’ (solid lines) 
given e = 1/3, the gas in Sagittarius B2 Main and North 


(dashed lines) and the Arches cluster from Espinoza, Selman 
fc Melnick] (|2009 1 (black, dash/dot) and Harfst, Portegies 


Zwart fc Stolte|(|2010[) (red, dash/dot). The plot shows that 


GO.253-1-0.016, ‘d’, ‘e’ and ‘f’ all contain enough mass such 
that they can form a YMG of M >10'* Mq given a SFE 
of 1/3. However, it is clear in all cases that the distribu¬ 
tion of mass is much less centrally concentrated than in the 
Arches, leading us to conclude that these clouds cannot form 
an Arches-like YMG at their current densities. Unless they 
were to somehow condense rapidly on global scales within a 
free-fall time (~ 0.5 Myr) before the onset of any widespread 
star formation, it seems implausible that these clouds will 
form clusters in accordance with scenario (i). We see that 
this is also true of the gas content of Sgr B2 Main and North. 
Despite the high mass and density of these regions, the gas 
is too extended on global scales. 

Figure 7 shows the resultant mass surface density pro- 
Hles as a function of radius for GO.253-1-0.016, ‘d’, ‘e’ and 
‘f’ given e = 1/3 (solid lines), the gas in Sagittarius B2 
Main and North (dashed lines), the proto-cluster(s) embed¬ 
ded within Sagittarius B2 Main and North as calculated 
from the UGHH region distribution (blue and red dash/dot, 
respectively) and the Arches cluster overlaid with the ht 
from [Espinoza, Selman fc Melnick (20091 (black open cir¬ 
cles, dash/dot). We again see that GO.253-1-0.016, ‘d’, ‘e’ 
and ‘f’ are less centrally concentrated than the Arches’ stel¬ 
lar distribution - further dynamical evolution would be re¬ 
quired to condense any resultant stellar populations if they 
were to form at the current densities of the clouds. We see 
that this is also true for the gas content of Sgr B2 Main 
and North. Even though they are significantly more mas¬ 
sive, dense and evolved than the other clouds, the gas is still 
smoothly distributed on global scales. This is in contrast to 
the inferred stellar content in Sgr B2, where the HII regions 
of both Sgr B2 Main and North are concentrated in a small 
volume (< 0.1 pc). This shows that dense, centrally concen¬ 
trated proto-clusters are able to form in clouds that are not 
very centrally concentrated, as predicted by hydrodynam- 


ical cluster formation simulations (Kruijssen et al. 20121. 


Furthermore, we find that, at least for this sample of clouds 
towards the Galactic centre, the hnal stellar distribution of 
an Arches-like YMG is more compact than the global distri¬ 


bution of stars and gas at any prior stage of the formation 
process. Given that these are the most massive and dense 
quiescent clouds yet found in the Galaxy, we therefore con¬ 
clude that scenario (i) is disfavoured as a likely mode of 
YMC formation. 

Having investigated the absolute mass and surface den¬ 
sity profiles, we now wish to identify the shape of the radial 
mass distribution so as to investigate any differences. Figure 
8 displays the normalised mass surface density profiles for 
GO.253-1-0.016, ‘d’, ‘e’ and ‘f’ (solid lines), the gas in Sagit¬ 
tarius B2 Main and North (dashed lines), the Sagittarius B2 
Main and North proto-clusters (blue and red dash/dot, re¬ 
spectively) and the Arches cluster (black, dash/dot), where 
all profiles have been normalised to unity at a radius of ~ 
0.3 pc (the resolution of the BGPS data at a distance of 8.4 
kpc). This plot shows that in all cases, the gas in the YMC 
progenitors is over-dense at large radii and under-dense at 
small radii, compared to the stellar components in both the 
Arches cluster and the Sgr B2 proto-clusters. This suggests 
that if these clouds are to form YMCs, the resultant stel¬ 
lar population would have to dynamically interact such that 
it would relax into a much more centrally-condensed dis¬ 
tribution. It is interesting that the gas distribution in both 
Sagittarius B2 Main and North on global scales is very sim¬ 
ilar to that in the quiescent ‘dust-ridge’ clouds. We see evi¬ 
dence for clustered massive star formation in Sgr B2, yet the 
global gas content looks identically distributed to that in a 
quiescent cloud. This is in conflict with scenario (i), which 
requires the gas to be more centrally concentrated prior to 
the formation of a YMC. 

It is also interesting to note the high mass surface den¬ 
sity of Sgr B2 proto-clusters on scales smaller than ~ 0.1 
pc. The bulk of the stellar mass, though roughly an or¬ 
der of magnitude lower than the Arches, is highly concen¬ 
trated within a small core region. It is difficult to determine 
whether this is consistent with scenario (i), in which case 
the stars may have formed in a centrally-condensed distri¬ 
bution, or whether it is consistent with scenario (ii), where 
the stars may be distributed in this way as a result of rapid 
dynamical interaction. 


4 CONCLUSIONS 


Using HiGAL far-IR continuum data and MALT90 mil¬ 
limetre spectral line data, we derive the global properties of 
four likely YMC precursors towards the Galactic centre. We 
hnd that these clouds, GO.253+0.016, ‘d’, ‘e’ and ‘f’, are all 
high mass (M ~ 10® Mq), cold {Tdust ~ 20 K) and dense (p 
~ 10* cm“®). They are all close to virial equilibrium and are 
therefore likely to be gravitationally bound. These results 
confirm that they are excellent YMC progenitor candidates 
- they are amongst the most massive and dense molecular 
clouds known to exist in the Galaxy, yet they are devoid of 
any widespread star formation. 

Comparison of the mass surface density distributions of 
these clouds with the Sagittarius B2 proto-YMC and the 
Arches YMC shows that these clouds are not compact or 
centrally concentrated enough to form an Arches-like YMC 
in their current state. If they are to form YMCs, dynamical 
evolution during the early formation must further condense 
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0.6 

0.44 
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d 
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8.4 
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19-23 

16.3 ± 1.5 
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1.4 

0.8 

0.3 

0.45 
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0.74 
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e 
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2.4 
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15.5 ± 0.5 

7.3 

0.7 

2.8 

0.9 

0.36 
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0.40 

14.5 

0.8 

f 

7.3 

8.4 

2.0 

18-22 

16.5 ± 3.2 

6.9 

0.9 

3.2 

0.8 

0.28 

4.3 

0.37 

15.4 
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Table 1. Global properties of clouds GO.253+0.016, ’d’, ’e’ and ’f’. The columns show mass (M), distance (D), radius (R), dust 
temperature (T^ust)? linewidth (AV), Virial mass Virial ratio (q) average volume density (n), average column density (Nh 2 )^ 

cloud-crossing time (tec), s ound-crossing time (tsc) -, free-fall time (t//), Mach number (Al) and the radius wit hin which 5 x 10^ Mq is 
enclosed (R5). * Result from Longmore et al. , 2012). ^ Galactrocentric distance estimate from 


Reid et al. 


to be at this distance. Sound speed is calculated using a gas temperature of 70 K \Ao et al.\201^ . 


'2009) - all clouds are assumed 



Figure 6. Enclosed mass as a function of radius for GO.253+0.016, ‘d’, ‘e’, ‘f’ given e = 1/3 (solid lines), the gas in Sagittarius B2 Main 
and North (dashed lines) and the Arches cluster from [Spinoza, Selman Sz Melnick] l |2009[ | (E09; black , dash/dot) and|Harfst, Portegies] 
|Zwart Stolte] ( |2010| | (HIO; red, dash/dot). GO.253+0.016* (red) ALMA cycle 0 + single-dish data ( [Rathborne et al.||2014b[ |. We see 
that in all cases, the gas in the proto-cluster clouds has a more extended distribution of mass, whereas the stellar population of the 
Arches is more centrally concentrated. 
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Figure 7. Mass surface density profiles as a function of radius for GO.253+0.016, ‘d’, ‘e’, ‘f’ given e = 1/3 (solid lines), the gas in 
Sagittarius B2 Main and North (dashed lines), the proto-cluster(s) embedded within Sagittarius B2 Main and North as calculated from 
the UCHII region distribution (blue and red dash/dot, respectively) and Arches cluster overlaid with fit from |Espinoza, Selman &: Melni^ 
l |2009[ | (E09; black open circles, dash/dot). These profiles show that the final stellar distribution of the YMC is more compact than the 
global distribution of stars and gas at any prior stage of the formation process. Note that the variation of the profiles for the HII region 
distribution in Sgr B2 at small radii is a result of the small number of sources detected there. The bump at R ~ 1.0 pc is due to Sgr B2 
Main entering the aperture at that radius. GO.253+0.016* (red) ALMA cycle 0 + single-dish data ||Rathborne et al.|2014b1l. 


the resultant stellar population. Furthermore, we find that 
the stellar content of both Sagittarius B2 Main and North 
is significantly more centrally concentrated than the global 
gas content. This shows that dense, centrally concentrated 
stellar clusters can form from gas that is not very centrally 
condensed, thus disfavouring a ‘monolithic collapse’ mode of 
YMC formation in which gas expulsion causes the YMC to 
end up less compact than the progenitor gas cloud. 
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Figure 8. Normalised mass surface density profiles as a function of radius for the gas in GO.253+0.016, ‘d’, ‘e’, ‘f’ (solid lines), the gas in 
Sagittarius B2 Main and North (dashed lines), the proto-cluster(s) embedded within Sagittarius B2 Main and North as calculated from 
the UCHII region distribution (blue and red dash/dot, respectively) and the Arches cluster from [fepinoza, Selman fc Melnick] | |2009[ l 
(E09; black dash/dot). Each profile has been normalised to unity at a radius of ~ 0.3 pc, as this is the resolution of the BGPS data. 
GO.253+0.016* (red) ALMA cycle 0 + single-dish data | |Rathborne et al.|2014bl |. The shapes of the profiles at larger radii suggests an 
evolutionary progression on global scales from a dispersed natal gas cloud to a centrally condensed stellar cluster as a function of star 
formation activity. 
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